Homogeneous and inhomogeneous magnetic phases of constrained dipolar bosons 
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We study the emergence of several magnetic phases in dipolar bosonic gases subject to three- 
body loss mechanism employing numerical simulations based on the density matrix renormalization 
group(DMRG) algorithm. After mapping the original Hamiltonian in spin language, we find a strong 
parallelism between the bosonic theory and the spin-1 Heisenberg model with single ion anisotropy 
and long-range interactions. A rich phase diagram, including ferromagnetic, antiferromagnetic and 
non-local ordered phases, emerges in the one-dimensional case, and is preserved even in presence of 
a trapping potential. 

PACS numbers: 67.85.Hj, 05.10.Cc, 75.10.Pq 



I. INTRODUCTION 

Recent experimental advances in controlling ultracold 
gases of magnetic atoms 1 - and polar molecules^ have 
paved the way to the investigation of several quantum 
many-body phenomenal. These setups naturally provide 
anisotropic, long-range dipolar interactions, which can 
be tuned and manipulated with high accuracy in order 
to access the physics of spin systems 4 - and Hubbard-like 
models^ loaded into optical lattices 5 -. Considerable theo- 
retical efforts have focused on one dimensional geometry, 
where non-local interactions play a fundamental role in 
stabilizing interesting phenomena such as supcrsolidity^, 
checkboard insulatori^ and insulating phases character- 
ized by non-local order parameter o 9 ' 10 . 

Furthermore, dissipative processes have emerged as 
an additional, relevant source of interaction. Two- 
body losses have been successfully employed to engineer 
hard core interactions in molecular gases 1 - 1 -, thus lead- 
ing to the stabilization of a Tonks-Giradeau gas, and 
three-body losses have been proposed as a dynamical 
source of three-body interaction^ 2 - which stabilizes dimer- 
supcrfluidit y 12 ! 13 , color-superfiuidity-i 4 - and Pfaffian-like 
state o 15 ' 16 with ultracold atoms. 

The aim of this paper is to investigate the interplay 
between local (two- and three-body) and non-local inter- 
actions in low dimensional systems of ultracold dipolar 
bosons. We focus on a one dimensional geometry, and 
find that the phase diagram of such systems strongly re- 
sembles that of the spin-1 Heisenberg model with Ising- 
like and single ion anisotropy-ii~— , or A — D model, 
extensively studied in the past in the contest of one- 
dimensional spin chains. We present numerical results on 
the phase diagram, which exhibits ferromagnetic, antifer- 
romagnetic and hidden order phases, and finally discuss 
the stability of these phases in presence of a trapping po- 
tential and density fluctuations, as naturally present in 
cold atomic and molecular setups. Several of the mag- 
netic phases discussed do not require strong dipolar in- 
teractions, and can thus be observed even with magnetic 



atoms, where dipolar interaction is usually much smaller 
than any other relevant energy scaled. 

The paper is organized as follows: in Sec. [IIJ we de- 
scribe the parallelism between constrained bosonic gases 
and spin systems and introduce the Hamiltonian, which 
is then investigated by DMRG simulations and strong 
coupling arguments in Sec. IIIIL In Sec. HVI we extend the 
numerical simulations to the inhomogeneous case; finally, 
we draw our conclusions in Sec. |V] 



II. BOSONIC HAMILTONIAN AND A-D MODEL 

Dipolar bosons confined in a one dimensional geom- 
etry and subject to a deep optical lattice are generally 
described by the following Hamiltonian^: 

H = -^&t& J - + -^n i (n i -l) + j u^n i + 

(id) i « 

% i<j u > 

Here, b\,bi,ni are bosonic creation, annihilation and 
number operator at site i, the first line describes the stan- 
dard Bosc-Hubbard model, where t is the hopping term 
between nearest neighbor sites and U the onsite interac- 
tion, and the last line includes trapping and long-distance 
dipolar potentials. The hopping coefficient t varies with 
the depth of the underlying optical lattice, whereas A can 
be tuned by varying the applied EC electric field; finally, 
the onsite interaction U depends on the short-distance 
details of the interparticle interaction^ 1 , and, for magnetic 
atoms, can be controlled by using Feshbach resonances^. 
The phase diagram of Eq. (Q]) with U > has been in- 
vestigated in several regimes: at unitary filling, a new 
insulating phase characterized by hidden order has been 
predicted between a Mott insulator and a charge density 
wav o 9 ' 10 , whereas for densities n < 1 and strong repul- 
sive interaction [/ » ( a devil's staircase of insulating 
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phases appears as a function of the chemical potential^.. 
However, the attractive regime U < has so far been 
neglected. This is partially due to the fact that losses 
given by strong three-body recombination are enhanced 
in this regime, thus making a time-dependent descrip- 
tion of the system more suitable in order to take into ac- 
count dissipative effects 2 ^. The situation can be though 
strongly simplified when the decay rate 73^ associated 
with three-body loss processes is much larger than the 
typical tunneling rate, i.e. 73 3> t: in this regime, a 
mechanism analogous to the quantum Zeno effect gives 
rise to an effective strong three-body repulsion, which can 
be implemented in the Hamiltonian with the additional 
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condition (bj) 

The opportunity to engineer strong three-body repul- 
sion has then two striking effects: i) the system is in 
general stable regardless of the sign of the couplings U, A 
and ii) the onsite Hilbert space is reduced to |0), |1), |2), 
thus resembling a spin-1 system. This correspondence 
is further clarified after introducing spin-1 operators 
, S~ , Sf and performing the following mapping: 



= 1-Sf 



b\ = aSr +P{SIS- +S-S1) (2) 



where a = (2+\/2)/4, /3 = -(2-^/4 are fixed by com- 
mutation relations, as described in appendix [A} From 
now on, we will consider a fixed density n = 1, then 
obtaining (fixing t = 1): 

U = - SfSr(J + JiSf + J 2 S] + J 3 S*S]) + (3) 

<i,j> 

sis* u 

(J - if + 2 



X0S?) 2 -fc£(i-£/2) 2 ^ 



where the first line includes a nearest-neighbor exchange 
with J — a 2 — /3 2 and correlated exchange terms with 
J\ = J2 = 2(3 and J3 = 4/3 2 , which break particle- 

hole symmetry, as required for constrained bosons. 

We notice that Eq. §5§ is a generalization of the so 
called A — D modet^, extensively studied over the last 
two decades both from analytical and numerical points 
of view. 

The spin-f A — D model presents a rich phase dia- 
gram: in addition to ferromagnetic and antifcrromagnetic 
(AFM) phases, in a broad region of the parameter space 
competition between local and non-local interactions fa- 
vors the so-called Haldane phase (as expected for inte- 
ger spin chains^), which displays a gap in the energy 
spectrum, a unique ground state (at least in the ther- 
modynamic limit, whereas it is four- fold degenerate for 
chains of finite size), a finite correlation length, and thus 
no long-range order even if it is possible to define suit- 
able string correlation functions that measure a hidden 
topological order. The spin liquid picture introduced by 
Tasakii^ provides a intuitive understanding of the Hal- 
dane phase: let us assign the presence of an effective spin- 
1/2 particle with spin pointing up (down) if at the i-th 
lattice site Sf = +!(—!) and no particles if 5f = 0. The 
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FIG. 1: (color online): homogeneous phase diagram for dipo- 
lar bosons on an optical lattice with three-body hard core 
constraint at filling n = 1 (see text): triangles, squares, dia- 
monds, black and red points denote numerical results, while 
the black dashed line describes an approximate strong cou- 
pling description for \U\ S> J. 



Haldane phase is then interpreted as a liquid in which 
these effective particles carry no positional order along 
the chain but still retain antifcrromagnetic (AFM) order 
in their effective spins. The positional disorder is asso- 
ciated with the absence of long-range order in the usual 
spin-1 correlation functions 
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x,y,z 



(4) 



whereas the spin- 1/2 magnetic order that we would get 
if all the sites with 5,f = were taken off from the chain 
is measured by the asymptotic value of the string order 
parameters (SOP)^: 



Oa(j) = (Si 
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x,y,z. 



(5) 



As shown thoroughly by Kennedy and Tasakii^ the A — D 
model possesses an hidden (non-local) Z2 x Z2 symmetry, 
and the non-vanishing values of the SOP can be under- 
stood as the breaking of such a symmetry. 



III. HOMOGENEOUS PHASE DIAGRAM. 

In order to exploit a complete parallelism between Eq. 
([3]) and the A — D model, we investigate its phase diagram 
in the homogeneous case, k = 0, by means of numerical 
simulations based on the density-matrix renormalization 
group (DMRG) algorithm 2 ^, truncating the dipolar in- 
teraction up to fifth-nearest-neighbors2£. Let us sum- 
marize the main results, as schematically presented in 
Fig. [U the A > region displays i) an antifcrromagnetic 
Neel-like phase (NP), where doubly occupied sites alter- 
nate with empty ones in a periodic pattern, ii) an Hal- 
dane insulator phase (HI), where doubly occupied and 
empty sites are separated by strings of single occupied 
one o 9 ' 27 (see Fig. [2]), and Hi) two superfluid phases, in 
which the superfluid components are single bosons (SF) 
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FIG. 2: (color online): cartoon of magnetic phases related to 
the model in Eq. (|ll3[) : from top to bottom, ferromagnetic, 
Haldane and Neel phase in bosonic and spin language(see 
text). 



and dimcrs (DSF) respectively. In the A < regime, 
both supcrfiuid phases collapse beyond a critical value 
of A into a ferromagnetic phase (FP), where the mu- 
tual attraction between bosons gives rise to a region of 
constant density n = 2. The system thus displays all 
phases and phase transitions of the A — D model with 
attractive single-ion anisotropy; there are however some 
quantitative differences. First, both SF and DSF, which 
correspond to the XY phases in spin language, extend 
on a broad region around A = due to the presence of 
correlated hopping terms which disadvantage long-range 
order. In addition, the HI region is present even at larger 
A, as expected due to long-range frustration of dipolar 
interactions with respect to antifcrromagnetic ordering^. 

Different phases are uniquely characterized by the 
asymptotic decay of correlation functions^. In the SF 
phase, both single particle and dimer superfluid correla- 
tions 

B(j) = (b\b i+j ) (x C x (j), V(j) = ((bj) 2 (b l+J ) 2 ) (6) 

decay algebraically; by contrast, in all other phases B is 
exponentially suppressed, whereas V decays algebraically 
in the DSF phase, as can be seen in Fig. [3] Mag- 
netic phases are instead characterized by a non- vanishing 
asymptotic value of certain correlation functions: in the 
HI, both O x , O z decay to a constant at long distances 
while C z vanishes exponentially, whereas in the NP O x 
decays exponentially and O z ,C z are constant. All mag- 
netic order parameters decay at long distances in both SF 
and DSF phases. Correlation functions have been com- 
puted by analyzing systems of size L=60, 80, 100 and 120 
sites, with up to 600 states per block, 4 sweeps and open 
boundary conditions. Fig. [4] describes typical decays in 
the SF (black, dashed), HI (red, thick) and NP (green, 
dot-dashed) of the magnetic order parameters. 

Haldane insulator - Neel phase . Any one of the pairs 
{C z ,O z }, {C z ,O x } or {O x ,O z } can be used to give an 
accurate description of the HI-NP boundary: from Figs. 
El El it can be inferred that the bulk asymptotic behavior 
of these correlators is well described already for L = 60. 
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FIG. 3: (color online): superfluid correlations in double log- 
arithmic scale as a function of the distance from the middle 
in a L — 120 chain. Red (dashed) and black (thick) lines rep- 
resent DSF and SF phase respectively, with U = —6.5, A = 
and U = -3, A = 0.05. 



This feature is not surprising, considering that in the A — 
D model this transition is believed to belong to the Ising- 
type universality class^. However, the same cannot be 
said when considering the other phase transitions present 
in the model. 

Superfluid - Dimer superfluid . The SF-DSF phase 
transition corresponds to a level crossing in the spec- 
trum between excitations with S z = ±1 and S z = ±2, 
with finite size gaps Ai and A2 respectively^. This con- 
dition stems from the fact that, beyond a critical attrac- 
tion — U C (A) 3> J, a finite energy is required to break 
dimer s^ and thus exciting the system in the S z = ±1 
sector would become energetically unfavorable. We de- 
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FIG. 4: (color online): magnetic order parameters in dou- 
ble logarithmic scale as a function of the distance from 
the middle in a L — 120 chain. Black (dashed), red 
(thick) and green (dot-dashed) represent (U = —3, A = 
0.05), (-3,0.7), (-3, 1.3) respectively. 
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FIG. 7: (color online). Left panel: dependence of A = Ai— A2 
on U for A = —0.1 and different chain lengths: L — 8 (black 
circles), L = 12 (red squares), L — 16 (green stars), L = 
20 (blue diamond) and L — 24 (orange triangles). Dashed 
lines are guides for the eye. Right panel: critical value of 
the SF-DSF transition as a function of 1/L for A = —0.1; 
circles represent numerical datas, line is a best fit of the type 
a% + a 2 /L + az/L 2 . 



FIG. 5: (color online). C z (x) correlation function at 
the HI-NP boundary for chains of different lengths L = 
60,80,100,120, panels a,b,c and d respectively; correlations 
are taken with respect to the center of the chain. Here, 
U — —3, and, from top to bottom, A = 1.1 (orange, dot- 
ted), A = 1.05 (blue, dashed), A = 1 (green, dot-dashed), 
A = 0.95 (red, thick) and A = 0.9 (black, dot-dot-dashed). 
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FIG. 6: (color online). String correlator O x (x) at the HI-NP 
boundary for chains of different lengths L — 60, 80, 100, 120, 
panels a,b,c and d respectively; correlations are taken with re- 
spect to the center of the chain. Here, U = —3, and, from top 
to bottom, A = 0.9 (black, dot-dot-dashed), A = 0.95 (red, 
thick), A = 1 (green, dot-dashed), A = 1.05 (blue, dashed) 
and A = 1.1 (orange, dotted). 



for periodic chains of several lengths and then by im- 
posing that, at the phase transition, the condition 
limt_ voo (Ai(i) — A2(L)) = is satisfied; a typical set 
of data for A = —0.1 is presented in Fig. [7j 

Superfluid - Haldane insulator and Dimer superfluid - 
Neel phase . Finite size calculations are also useful to 
better shape the SF-HI and DSF-NP transitions, which, 
in analogy with the A — D model^^r— , should belong 
to the Bcrezinskii-Kostcrlitz-Thouless (BKT) universal- 
ity classic, albeit the non trivial nature of non-local inter- 
action can in principle lead to different critical behaviors. 
BKT transitions are usually hard to determine due to the 
exponential opening of the gap; however, string order pa- 
rameters have been shown to provide a rather accurate 
estimate of the transition point a 29 ! 31 . In the following, 
we consider systems with periodic boundary conditions, 
in order to avoid boundary effects, with up to L = 42 
sites, and calculate the string order parameter from the 
first to the L/2 + 1 site, O z (x — x' — L/2), which rep- 
resents a suitable order parameter for both SF-HI and 
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tcrminc the phase boundary (marked by blue diamonds 
in Fig. [T]) by calculating finite size gaps 
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FIG. 8: Finite-size scaling of O z (L/2) for U = -0.8 by 
using Eq. ([7| (see text). From top to bottom: A = 
0.5, 0.45, 0.4, 0.35, 0.3, 0.25, 0.2, 0.15 and 0.1. 
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FIG. 9: (color online). Upper panel: string order parameter 
d for U = —0.8, —2 (black circles and red squares respec- 
tively) as a function of A. Lower panel: Ci for A = 0.1, 0.2, 0.3 
(black squares, red triangles and blue circles respectively) as 
a function of U. In both panels, lines are guide for the eye, 
and the size of each point denotes the maximum error on the 
extrapolated value Ci. 



DSF-NP transitions^. Then, we estimate its asymptotic 
value C\ by fitting the data with the following scaling 
form: 

0,(L/2) = C 1 + ^ r . (7) 

This method has been successfully employed to study the 
same transition in the \—D model within an exact diago- 
nalization approach^. The usual error of this procedure 
is related to the DMRG truncation error, always smaller 
than 3 * 10~ 5 , and to the algebraic fh— : by employing 
numerical datas with 20 < L < 42, we then estimate that 
the asymptotic value C\ is non-vanishing within numer- 
ical error as long as C\ > 0.005. A typical example of 
the estimate of C\ is described in Fig. [51 where differ- 
ent datas are presented for U = —0.8,0.1 < A < 0.5. 
In Fig. [9j we plot C\ as a function of both U (fix- 
ing A = 0.1,0.2,0.3) and A (fixing U = -0.8,-2); the 
asymptotic value of the string order parameter increases 
with both A and \U\. The corresponding transition points 
are marked by red triangles in Fig. [TJ due to the small 
system sizes analyzed here, the numerical data presented 
above can be considered as an approximate estimate on 



the BKT transition points, an accurate one needing more 
specific technique o 20 ' 31 . 

Ferromagnetic phase - Superfluid . A clear evi- 
dence of the FP-SF transition emerges instead when the 
ground state energy density £qs{L) approaches the exact 
value for ferromagnetic states £ fp in the thermodynamic 
limits. The transition line, marked by black dots in 
Fig-©7 is obtained by requiring that lim^^oo Egs{L) = 
£fp- 

A. Strong coupling regime 

In the large \U\ regime, the quantitative difference 
between dipolar and ncarest-neighbor(NN) interaction 
can be investigated with a perturbative argument. If 
\U\ ^s> |A|, J, the effective Hilbert space is reduced to 
S z = ±1, that is, bosons are tightly bound in dimers, 
so that we can map the spin-1 problem into a spin- 1/2 
theory employing the following identities^: 

Sf = 2 s z v SfSf = 2sj, Srsr = 2sj (8) 

S+S- = 2(1/2 + sf), S-S+ = 2(1/2 - sf ) (9) 

where s is a spin- 1/2 operator. After a proper reseat- 
ing, the strong coupling Hamiltonian including only NN 
interaction is mapped into a spin-1/2 XXZ chain: 

Hsc = £ ( S + S 7 + (1-A)< S |), A = -2/A\U\. (10) 

<i,j> 

From the exact solution of Eq. ([TUp i. we argue that the 
system is in a DSF phase as long as — < A < 0; 
the DSF-FP and DSF-NP transitions are located at 
A-NN = — WT ' respectively. We can now compare this 
criterion, derived considering only NN interactions, with 
the numerical one, obtained from DMRG as previously 
described. The DSF-FP transition (black dashed line in 
Fig HI ) is then predicted at —A ~ 2/\U\, whereas nu- 
merical values (including dipolar interaction) for, e.g., 
U = —10, indicate —A ~ 1.75/|E/|; in this regime, dipo- 
lar interactions show a small quantitative difference with 
respect to standard NN couplings. 

IV. EFFECT OF A TRAPPING POTENTIAL. 

The observation of the different magnetic orders dis- 
cussed above in a standard cold atom experiment is 
strictly related to the possibility of stabilizing these 
phases even in an inhomogeneous background. In fact, 
atoms and molecules are loaded into a trapping poten- 
tial, which introduces a position-dependent term in the 
Hamiltonian with a minimum at the trap center, namely, 
the first term in the second line of Eq.(TT]). In this setup, 
another energy scale comes into play; particles will try 
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to minimize their potential energy by concentrating in 
the middle of the trap, thus displaying a strong spatial 
dependence of the local density (rii). This feature is in 
sharp contrast to the ideal configuration needed to re- 
alize the magnetic phases described above, all of them 
requiring a constant density in the thermodynamic limit. 
Our goal here is to investigate what are the proper trap 
configurations needed to stabilize a magnetic phase in 
this inhomogeneous setup: in particular, we will focus 
our attention to the region in the middle of the confining 
potential, where it is usually easier to create large re- 
gions of space at constant density^^. First of all, let us 
briefly discuss what happens the in the FP: in analogy to 
a standard Mott-like phased, it can always be realized by 
considering a sufficiently strong trap such that the den- 
sity is maximized in the middle, (rii) = 2. This simple 
argument cannot be extended to neither HI or NP: in 
fact, a very strong trap will simply destroy these types of 
order. In this section, we will thus focus on the stability 
of these two orders in presence of a trapping potential. 

We identify a certain magnetic order in a region of 
space by requiring that i) the region is at constant den- 
sity, (rii) = 1, and ii) the proper order parameters with 
respect to the middle of the trap, defined as: 



O a=x , z (L/2,j) = (S2 /2 e l ^^< 3+ ^^S2 /2+J ) (11) 



si 



C z (L/2,3) = (-iy(Sl /2 Sl /2+J ) 



(12) 



where j is the distance from £/2, behave as expected in 
the HI or in the NP up to a certain range. We performed 
DMRG simulations on a L = 80 sites chain, fixing as 
energy unit for the trapping potential k* = l/(L/2) 2 = 
1/1600, and kept as much as 600 states per block with 10 
finite-size sweeps^. Since a trapping potential favors a 
configuration where the particles are in the middle of the 
chain, we focused on a N = 40 particles system: this as- 
sures that the density close to the chain boundary rapidly 
goes to zero, thus avoiding possible finite-size effects and, 
at the same time, allows for a constant density of order 
1 in the middle of the system. 



A. Haldane order 

We start our treatment by considering the possibility 
to stabilize hidden order in an inhomogeneous system. 
As a sample configuration, we fixed U = —2.5, A = 0.9 
such that the corresponding homogeneous phase at inte- 
ger filling is an HI. In a very shallow trap, the non-local 
interparticlc repulsion would drive the system in a very 
dilute limit with (rii) < 1 all over the trap, whereas in 
the opposite strong trap limit, an high density region 
with (rii) > 1 would be stabilized in the middle. We 
shall then focus on an intermediate regime in order to 
satisfy the density requirement (rii) ~ 1 needed in the 
HI. 




FIG. 10: (color online): density distribution as a function of 
the distance from the trap center for U = —2.5, A = 0.9. Top 
panel: fixed population N = 40 and different trap strength 
k — 7.5k* (black, dotted), 8.5fc* (red, thick), 9.5fc* (green, 
dashed) and 10.5A:* (blue, dot-dashed). Bottom panel: fixed 
trap strength k = 9k* and different populations: N=36 
(black, dot-dashed), 38 (red, dashed), 40 (green, thick), 42 
(blue, dot-dot-dashed) and 44 (orange, dotted). 
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FIG. 11: (color online): magnetic order parameters C z (L/2,j) 
(top panel) and O z (L/2, j) (bottom panel) as a function of the 
distance from the trap center for U — —2.5, A = 0.9, N = 40 
and different values of k: k — 7.5k* (black, dotted), 8.5fc* 
(red, thick), 9.5fc* (green, dashed) and 10.5fc* (blue, dot- 
dashed). 



In Fig. 1101 upper panel, we plot the density distribu- 
tion as a function of the distance from the trap center for 
different values of k, 7.5 < k/k* < 10.5: the requirement 
(rii) = 1 is satisfied for values of the trap strength inside 
the interval 8.5 < k/k* < 9.5. In order to verify whether 
hidden order is present or not, we plot the relevant mag- 
netic order parameters as defined in Eqs. 111! 1121 in Fig. 
QTJ in the interval 8.5 < k/k* < 9.5, the string order 
parameter O z is constant up to a certain distance from 
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FIG. 12: (color online) : magnetic order parameters C z (L/2,j) 
(top panel) and O z (L/2,j) (bottom panel) as a function of the 
distance from the trap center for U = —2.5, A = 0.9, k — 9k* 
and different number of particles N=36 (black, dot-dashed), 
38 (red, dashed), 40 (green, thick), 42 (blue, dot-dot-dashed) 
and 44 (orange, dotted). 



the trap, and at the same time C z decays, then proving 
that particles close to the trap center display HI; outside 
of the constant density region, the order is lost, as can 
also be seen by looking at O z . 

However, while an accurate fine tuning of the trap 
strength does not present major difficulties in a typical 
experimental setup, a proper control over populations in 
a tube is challenging, and it is thus instructive to investi- 
gate small population unbalance with respect to the pre- 
vious N = 40 case. In Fig. [101 lower panel, we plot the 
density distribution at a fixed trap strength k = 9k* for 
different total number of particles N = 36, 38, 40, 42 and 
44, while the corresponding order parameters are plotted 
in Fig. [12j We notice that the Haldane phase is unstable 
when N > 42 since too many particles concentrate in the 
middle of the trap, whereas it is stable for TV < 40; we 
can then conclude that a large population difference of 
order 5N ~ 0.1 prevents the HI phase to stabilize in the 
center of the trap. 



B. Antiferromagnetic ordering 
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FIG. 13: (color online): density distribution as a function of 
the distance from the trap center for U = —4, A = 1.5. Top 
panel: fixed population N = 40 and different trap strength 
k = 5k* (black, dotted), 7k* (red, thick), 9k* (green, dashed), 
Ilk* (blue, dot-dashed) and 20fc* (orange, dot-dot-dashed). 
Bottom panel: fixed trap strength k — 12k* and different 
populations: N=36 (black, dot-dashed), 38(red, dashed), 40 
(green, thick) and 42 (blue, dot-dot-dashed). 
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FIG. 14: (color online): magnetic order parameters C z (L/2,j) 
(top panel) and O z (L/2, j) (bottom panel) as a function of 
the distance from the trap center for U = —4, A = 1.5, N = 40 
and different values of k: k — 5k* (black, dotted), 7k* (red, 
thick), 9k* (green, dashed), Ilk* (blue, dot-dashed) and 20fc* 
(orange, dot-dot-dashed). 



We turn now our attention to the NP by fixing U = 
—4, A = 1.5. As already discussed for the HI, a very shal- 
low trap is not sufficient to stabilize antiferromagnetic 
order in the trap due to density requirements, whereas 
a too strong trap would prevent it by concentrating too 
many particles in the trapping potential minimum. The 
density distribution as a function of the trap strength in 
the interval 5 < k/k* < 20 is presented in Fig. [T31 up- 
per panel; a large region with (n.;) = 1 is stable in the 
middle as long as k > 5k*, and, remarkably, the size of 
this region increases with increasing k, including up to 
30 particles when k = 20k* . The corresponding magnetic 



order parameters are plotted in Fig. [T4l both C z and O z 
are constant in the middle of the trap as long as k > 5A;* , 
and their plateau extends all over the constant density 
region. 

Small changes in the total number of particles do 
not alter this picture significantly. In Fig. Q2J lower 
panel, we plot the density distribution at a fixed trap 
strength k = 12k* and different total number of particles 
N = 36, 38, 40, 42; a constant region in the middle of the 
trap is always present and, in addition, it displays an- 
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FIG. 15: (color online) : magnetic order parameters C z (L/2,j) 
(top panel) and O z (L/2, j) (bottom panel) as a function of 
the distance from the trap center for U = —4, A = 1.5, k — 
12k* and different number of particles: N=36 (black, dot- 
dashed), 38(red, dashed), 40 (green, thick) and 42 (blue, dot- 
dot-dashed). 



tifcrromagnctic correlations, as can be inferred from the 
magnetic order parameters presented in Fig. 1151 

We can then conclude that, while both the HI and NP 
are incompressible, the former requires a finer tuning of 
the trapping potential and a more accurate control of the 
population of the system in order to be stabilized in the 
centre of the trap. 



V. CONCLUSIONS 

We have investigated the effect of a three-body hard- 
core constraint in a one dimensional system of dipolar 
bosons such as magnetic atoms or polar molecules con- 
fined in a one dimensional tube by optical lattices. After 
mapping the original problem in spin language, a strong 
parallelism between the system and the A— D model is es- 
tablished and confirmed by DMRG calculations; we have 
shown that ferromagnetic, antifcrromagnctic and hidden 
orders can be stabilized in this framework, and that dipo- 
lar interactions give rise to small quantitative changes in 
the phase diagram with respect to more usual nearest- 
neighbor interactions. 

Finally, we have considered the effect of a confining 
potential, as usually present in cold atomic gas experi- 
mental setups. Remarkably, both HI and FP can be sta- 
bilized in a large region in the middle of the trap even if 
the total population is not exactly controlled; the realiza- 
tion of the HI needs however an appropriate trap tuning. 
This feature opens up the possibility to observe strongly 
correlated phases in cold gases of magnetic atoms such as 
Dy, Cr or £r— i 3 -, which are usually characterized by rela- 
tively small dipolar interactions. All of these phases can 
be probed via noise correlations^, or, in the HI case, via 
Bragg spectroscopy^ or in-situ imagin g 35 i 37 . Finally, this 



setup can be adapted to investigate spin-1 Hcisenberg- 
like models in 2-D systems, where various interesting 
phases such as field induced supersolidity have been re- 
cently suggested 3 ^ or else extended in order to consider 
the effect of disorder in such systems. 
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Appendix A: Constrained bosons - spin-1 mapping 

The constraint |0) = allows us to make the 

following correspondence between the reduced bosonic 
Hilbert space and that of a spin-1 

|o>-Ht) |i>-H6) |2>-H+) (ai) 

where | f), |0) and | I) are eigenstates of S z with eigen- 
values +1, and —1 respectively. The corresponding 
operator mapping is: 

tfb = 1 - S z (A2) 

b = aS + + p (S z S + + S + S z ) (A3) 

where the coefficients a, /3 have to be determined by im- 
posing the correct action on the Hilbert space and com- 
mutation relations. Verifying the former, we have: 

6|0) = [aS+ + f3(S z S+ + S+S z )]\t) =0 
fc|l) = |0) [aS + +f3(S z S+ + S+S Z )} |0) (A4) 
= (cn/2 + /V2)|t) 

6|2)=V2|1) [aS++l3(S z S+ + S+S z )]\i) 

= (cV2-/V2)|0>. (A5) 

It follows then 

a + fi = —= a -,3 = 1. (A6) 

Furthermore, if we write down the number operator in 
terms of spin-1 operators 

6*6 ->S = a 2 S~S + + f3 2 (S~S Z S + S Z + S~S Z S Z S++ 
+ S Z S-S + S Z + S Z S-S Z S + ) (A7) 
+ a(3(S-S + S z + 2S-S Z S + + S Z S-S + ) 

and we apply it to number eigenstates, we get: 

6+6|0> = S\ t) = 

&t&|l) = |l) S\0) =2(a + /3) 2 |0) = |0) (A8) 
6t&|2> = 2|2) S\i) =2{a-P) 2 \i) = 2| |). 
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We will now show that the operators defined by our map- 
ping satisfy the correct commutation relations. Since 
we are considering constrained bosons the usual bosonic 
commutation relation becomes^ 

[6,6+] =|0)(0| + |1)(1|-2|2}(2|. (A9) 

Recalling that 

[S + , S~] = 2S Z (A10) 
= 2 (s (S + 1) - (S z ) 2 ^j 

we have 

[b, 6*] = -8aj3 + 2 (a 2 + p 2 ) S z + 12a/3 {S z f (All) 



and then 

(0| [6,&t] |0)=2(a + /?) 2 = l 

(1| [b,tf] |1) = -8a/? = 1 (A12) 

(2| [6,fet] |2)=-2(a-/3) 2 = -2. 
In order these relations to be satisfied we must have 

2 + \/2 2-V2 
a= — - — /3 = — . (A13) 
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